Lagrangian time correlations of vorticity alignments in isotropic turbulence: 

observations and model predictions 
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Motivated by results from recent particle tracking experiments in turbulence (Xu et al., Nat. 
Phys. 7, 709 (2011)), we study the Lagrangian time correlations of vorticity alignments with the 
three eigenvectors of the deformation-rate tensor. We use data from direct numerical simulations 
(DNS), and explore the predictions of a Lagrangian model for the velocity gradient tensor. We find 
that the initial increase of correlation of vorticity direction with the most extensive eigen-direction 
observed by Xu et al. is reproduced accurately using the Lagrangian model, as well as the evolution 
of correlation with the other two eigendirections. Conversely, time correlations of vorticity direction 
with the eigen-frame of the pressure Hessian tensor show differences with the model. 
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In a recent communication [l| , the incrtial range vortic- 
ity vector defined as the rotational motion of an ensemble 
of four Lagrangian particles in a turbulent flow ( "tetrad" 
0) has been shown, at short times, to display growing 
alignment and correlation with the eigenvector associated 
to the most positive eigenvalue of the associated initial 
deformation tensor. It is well established that at 

any given time the alignment between vorticity is most 
likely to be with the intermediate strain-rate eigenvec- 
tor. However, how said correlation changes as function of 
time-delay between vorticity and strain-rate is an aspect 
of turbulence fine-scale structure that had not received 
much attention before. The time evolution of alignment 
correlations provides an interesting observable on which 
turbulence models can be tested. The evolution of align- 
ments between material lines and vorticity, and viscous 
mechanisms of tilting of vorticity, have been studied in 
experiments and simulations 0] • 

In this Letter, we examine the Lagrangian time cor- 
relation of vorticity and the various eigenvectors of the 
deformation rate tensor, as predicted by a Lagrangian 
stochastic model Q . Results are compared with the 
results of [l| and additional data from Direct Numeri- 
cal Simulation (DNS) of the Navier-Stokes equations in 
moderate Reynolds number isotropic turbulence. In the 
same spirit as the analysis proposed in Ref. we also 
study the short time alignment of vorticity with the pres- 
sure Hessian eigen- vector directions. The pressure Hes- 
sian tensor is a key quantity in the transport equation of 
the velocity gradients A, where Aij = dui/dxj and u is 
the velocity vector, given by 



dt 



/AA 
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where d/dt is the Lagrangian time derivative, Pij ~ 
dp/dxidxj the pressure Hessian (divided by fiuid den- 
sity) and ly the fiuid's kinematic viscosity. The results 
obtained from DNS are compared with thepredictions of 
the Lagrangian stochastic model of Ref. [Si . 

As in Ref. we focus now on the Lagrangian time 
correlation of vorticity a; = V A u with eigen-vectors 



of the rate of deformation tensor S = (A + A )/2. We 
define, as in Ref. the unit-norm vorticity vector direc- 
tion e^{t) = u>{t)/\uj{t) \ and the orthonormal eigenframe 
{ei{t), e2(t), e3(t)) of S associated to its 3 ordered eigen- 
values Xi{t) > X2{t) > A3(t). Incompressible fluids are 
considered, so that Xi{t) + X2{t) + X3{t) = which implies 
that Ai(i) > and A3(i) < 0. Following Ref. [H, who 
focused on the most extensive eigen-direction (ei(i) with 
1 = 1), the following Lagrangian correlation function is 
considered 



C^{T) = {[e,[t).e^{t + T)f) 



(2) 



for i € {1,2,3} and the temporal displacement r (and 
averaging in time t) is performed along Lagrangian tra- 
jectories. 

We consider predictions from the Lagrangian model 
developed in Ref. Q for the velocity gradient tensor 
in turbulent flows and given explicitly by the following 
stochastic differential equation 



dA = -A^ 



Tr(A^ 
Tr(C7 
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A\dt + dW 



(3) 

The second and third terms in the right hand side of 
Eq. 13] arc closures for, respectively, (minus) the pressure 
Hessian and viscous Laplacian that govern the time evo- 
lution of A, as it can be seen by comparing with Eq. [1] 
(see Ref. for a review of this and other models) . The 
term W is a tensorial delta-correlated noise term that has 
been added in order to represent possible forcing effects, 
e.g. from neighboring eddies. The "recent Cauchy-Grcen 
tensor" , which arises after invoking the "recent fluid 
deformation" approximation Q , can be expressed in 
terms of matrix exponentials: 

= e^'-^e^"^'" , (4) 

where is the Kolmogorov time-scale (see Ref. for 
details). Hence, the modeled pressure Hessian tensor 
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-C.,. can be highly anisotropic, mirroring 
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FIG. 1: (a) and (e) Correlation Ci{T) (Eq. 0) of the vorticity direction etj(t) and the eigenframe of the deformation rate, ei{t) 
in DNS flows and from the model (Eq. |3]): dashed line for maximal straining direction (i.e. i — 1), solid line for intermediate 
straining direction (i = 2), and dot-dashed line for most contractive eigendirection {i — 3). (b,c,d) and (f,g,h) Probability 
Density Functions (PDFs) of \ei{t).euj(t + t) \ at various time lag r represented by different symbols, i.e. = To < ttd < To < rv 
and indicated in (a) and (b). 



the deformation undergone by the fluid as represented 
by the inverse of Cr^. ■ If one is interested in the velocity 
gradient tensor of a 'coarse-grained' velocity field (such 
as occurs in large-eddy simulations at length-scale A, or 
in the context of "tetrads" when they span a typical scale 
A in the inertial range) , the Lagrangian stochastic model 
may also be interpreted as a model for the coarse-grained 
velocity gradient tensor if the Kolmogorov time-scale 
is replaced by the corresponding eddy turn-over time at 
the appropriate scale ta = g-i/3/\2/3^ where e is the 
mean dissipation. In Eq. [31 the viscous term includes 
the time-scale T, corresponding to the Lagrangian inte- 
gral time-scale of the velocity [8|. 

The model (Eq. consisting of 9 (8 independent) 
stochastic differential equations, can be run with arbi- 
trary initial conditions and it generates stationary statis- 
tics for the velocity gradient tensor elements. In partic- 
ular, signals corresponding to the time histories of the 
antisymmetric part of Aij (the vorticity vector and its 
direction Buit)) and of the directions of the strain-rate 
eigenvectors ei{t) are readily obtained from the model 
runs. The correlation functions are then evaluated by 
averaging signals over long (10^ T) records of model sig- 
nals. 

In the following, as comparison we will also make 
use of data from a standard direct numerical simulation 
(DNS) of the Navier-Stokes equations. DNS is based 
on a pseudo-spectral (de-aliased according to the |- 
rule) method with 2nd-order accurate Adams-Bashforth 



time stepping; the computational box is cubic (size 2tt) 
with periodic boundary conditions in the three directions 
and spatial resolution of 256'^. Statistical stationarity is 
maintained by an isotropic external force acting at low 
wavenumbers in order to ensure a constant power injec- 
tion. It provides, in the units of the simulation, a con- 
stant energy injection rate of e = 0.001. The kinematic 
viscosity of the fluid isiy ~ 0.0004. The Kolmogorov scale 
is rjK = 0.016 so that dx/rjK ~ 1-5 (with dx — 27r/256). 
The Taylor-based Reynolds number is of order R\ ^ 125. 

The Lagrangian stochastic model is run with the pa- 
rameter Tk/T — 0.1, which is appropriate for modeling 
turbulence with Rx ^ 150 We compare in Fig. [TJa) 
and (e) the correlation functions Ci (r) as predicted from 
the model runs (Eq. [3]) and measured from DNS. For 
both the DNS and the model, time is normalized by the 
"integral" correlation time of a single tensor compo- 
nent (e.g. ^11 ), namely 

1 /■+°° 

r,, = (All (t) An {t + T))dT . 

It is found that in DNS, ~ 2rfc with = \fvje and 
for the model, « 1.2Tfc with r^/T = 0.1. We see first 
that the general trends observed in DNS, and as reported 
in Ref. [l| from DNS and for inertial-range tetrads, are 
reproduced quite well by the stochastic model (Eq. [3]). 
More specifically, for time lags t < r,, we see that ei and 
e.^ are increasingly correlated, meaning that we observe. 
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at short times, a more and more pronounced alignment of 
vorticity with the most extensive eigen-direction of the 
strain-rate tensor. This was one of the main observa- 
tions reported in Ref. The general trend of initial 
alignment with the most extensive strain-rate direction 
is expected from vortex stretching, as discussed (e.g.) in 
Majda (1991) [l^l- Conversely, C3(t) decreases at short 
times, i.e. vorticity decorrelates with the most contrac- 
tive direction. As far as the intermediate eigen-direction 
is concerned, C2(t) decreases monotonically towards the 
uncorrelated limiting value Ci{T) = 1/3. The model 
shows a slightly faster decay compared to the DNS, in 
the units of r^. 

As proposed in Ref. [l|, we also display in Figs. 
[Ii;b,c,d) for DNS, and in Figs. [Ii;f,g,h) for the model, 
the probability density functions (PDFs) of the cosine of 
the angles between the direction of vorticity at time t + r 
and the three initial eigen-directions of the deformation 
at time t. Remarkably, the trends and results observed 
in DNS are reproduced quite well by the model. At van- 
ishing time-lags r = (symbols o), is uncorrelated 
with ei, mostly aligned with 62 and mostly orthogonal 
to 63. Then, for the first non vanishing time lag tq, we 
observe that e^{t + Tc\) becomes more aligned with ei(t), 
less aligned with e2(<) and more orthogonal to 63 (t). For 
later time lags < ry , alignments of vorticity with the 
eigenframe relax towards the uncorrelated situation, in 
which the PDF is flat. 

We see thus that the alignments of vorticity with the 
eigen-frame of the deformation rate tensor are well pre- 
dicted by the stochastic Lagrangian model (Eq. |3|). A 
key reason is that the model contains the exact "veloc- 
ity gradient self-stretching term" — in its right hand 
side. This term includes the familiar vortex stretching 
mechanism that determines the antisymmetric part of 
the tensor and hence the evolution of vorticity. There- 
fore, it may be expected that other Lagrangian stochastic 
models such as the model with prescribed log-normal dis- 
sipation [nj , the "tetrad model" 0] or the "Lagrangian 
linear diffusion model" [l^l that include the velocity gra- 
dient self-stretching term, should also display the trend 
of growing alignment of vorticity direction with the most 
extensive strain eigen-direction. 

Another important and related question is the time 
correlation of vorticity with the pressure Hessian eigen- 
frame. The pressure Hessian P is a key quantity that 
enters into the Lagrangian dynamics of the velocity gra- 
dient tensor A (Eq. [T|) [5| . As underlined in the literature 
understanding the time evolution of the align- 
ments of vorticity with the eigenframe of the deviatoric 
part of the pressure Hessian P'' = P — tr(P)I/3, in both 
the Euler and Navier-Stokes equations, is important. It 
has been reported that in stationary turbulent flows Q , 
vorticity gets preferentially aligned with the eigenvector 
associated to the intermediate eigenvalue of P''. This 
property was found well reproduced by the stochastic 
model. 

We display in Fig. [5] the Lagrangian time behavior of 



alignments of vorticity with the eigen-frame of P , in a 
similar way as in Fig. [1] but now with the deformation- 
rate tensor S replaced by P''. We define the correlation 
function 

A(r) = ([b,(t).e„(t + T)]') , (5) 

where hi(t) for i G {1, 2, 3} are the eigenvectors of P'' and 
the time is, again, understood to be along Lagrangian 
trajectories. In DNS (Fig. [H^a)), we see that at short 
time (i.e. r < r^), Di (resp. D3) decreases (resp. in- 
creases), in the opposite way as it was observed for the 
deformation-rate tensor (the C^'s). As far as the inter- 
mediate pressure Hessian eigendirection is concerned, D2 
increases slightly. For larger time lags, i.e. t > Ttj, all 
the functions -Di(T) relax towards the uncorrelated value 
Diij) = 1/3. In the stochastic model (i.e. in Fig. [DJe)), 
it seems that even at short times, the functions Diir) 
tend immediately from their initial values towards the 
uncorrelated value of 1/3: there is no intermediate non- 
trivial short time behavior. 

In Figs. [2l^b,c,d) (resp. Figs. [2][f,g,h)), we display the 
PDFs of the angle of vorticity with the hi for various time 
lags. At vanishing time lag Tq = 0, the PDFs are similar 
to the ones already displayed in Ref. In particular, in 
both DNS and in the model, the preferential alignment 
of vorticity and b2 is confirmed. See Ref. Q for further 
interpretations and discussions about this behavior. As 
seen in Figs. [2ljb,d) and Figs. [2l^f,h), at short times, 
vorticity becomes orthogonal to bi and becomes aligned 
to b3. This is not reproduced by the model in which 
we observe a simple relaxation towards the uncorrelated 
case (all the PDFs become flat). Recall that the model is 
based on assuming that the upstream Lagrangian pres- 
sure Hessian is isotropic, and that the tensor A does not 
change during the (recent) time tk- These assumptions 
are likely to be violated in real turbulence thus causing 
the observed differences. Note, however, that the evolu- 
tion of vorticity alignment does not directly depend upon 
the (symmetric) pressure Hessian tensor, but only indi- 
rectly through its dependence on strain-rate. Thus, elu- 
cidating the real causes of the observed model limitations 
is challenging and beyond the scope of this Letter. 

To conclude, we have shown in this Letter that at short 
times, vorticity becomes aligned with the most extensive 
eigendirection. This agrees with the trends observed for 
tetrads obtained at larger scales in the experimental and 
DNS investigation in Ref. [l| • This behavior is quite con- 
sistent with linear stretching of the vorticity field. The 
observed behaviors can be successfully predicted by the 
Lagrangian stochastic model based on a closure of the 
Pressure Hessian using the Recent Fluid Deformation 
approximation We have furthermore extended the 
study to the alignments of vorticity with the eigen-frame 
of the pressure Hessian, and we have underlined what is 
and what is not predicted by the model. 

We are very thankful to Emmanuel Leveque for provid- 
ing us with the DNS data. Computations have been per- 
formed by using the local computing facilities (PSMN) at 
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FIG. 2: (a) and (e) Intercorrelation Ci(r) (Eq. [2]) of the vorticity direction e„(f) and the eigenframe hi{t) of the deviatoric 
part of the pressure Hessian P'' in DNS fiows and from the model (Eq. [3]): dashed Une for maximal eigenvalue /3i, solid line 
for intermediate eigenvalue 132), and dot-dashed line for most negative eigenvalue Ps- Panels (b,c,d) and (f,g,h) are PDFs of 
\hi{t).eu:{t + t)\ at various time lag r represented by different symbols, i.e. = ro < thi < To < tv and indicated in (a) and (b). 
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